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QCD sum rules are analyzed with the help of the Maximum Entropy Method. We 
develop a new technique based on the Bayesion inference theory, which allows us to directly 
obtain the spectral function of a given correlator from the results of the operator product 
expansion given in the deep euclidean 4-momentum region. The most important advantage 
of this approach is that one does not have to make any a priori assumptions about the 
functional form of the spectral function, such as the "pole + continuum" ansatz that has 
been widely used in QCD sum rule studies, but only needs to specify the asymptotic values of 
the spectral function at high and low energies as an input. As a first test of the applicability 
^ \ of this method, we have analyzed the sum rules of the p-meson, a case where the sum 

C^ , rules are known to work well. Our results show a clear peak structure in the region of 

the experimental mass of the p-meson. We thus demonstrate that the Maximum Entropy 
Method is successfully applied and that it is an efficient tool in the analysis of QCD sum 
rules. 
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§1. Introduction 



The technique of QCD sum rules is well known for its ability to reproduce various 
properties of hadrons.'2''l2li Using dispersion relations, this method connects pertur- 

j-C^. bative and nonperturbative sectors of QCD, and therefore, allows one to describe 

inherently nonperturbative objects such as hadrons by the operator product expan- 

CN ' sion (OPE), which is essentially a perturbative procedure. The higher-order terms of 

the OPE contain condensates of various operators, which incorporate information on 

V/^ \ the QCD vacuum. Hence, QCD sum rules also provide us with nontrivial relations 

3^ ' between the properties of hadrons and the QCD vacuum. 

Since the early days of the development of QCD sum rules, the range of appli- 

^^ \ cations of this method has been constantly expanding, which has helped to explain 

f^ ' many aspects of the behavior of hadrons. Nevertheless, QCD sum rules have always 

been subject to (justified) criticism. One part of this criticism is of mainly technical 
nature, pointing out that the analysis of QCD sum rules often is not done with the 
necessary rigor, namely, that the OPE convergence and/or the pole dominance con- 
dition are not properly taken into account. Many of the recent works that followed 
C^ \ the claimed discovery of the pentaquark 6*^(1540) are examples of such a lack of 

rigor. Nonetheless, these technical problems can be overcome if the analysis is done 
carefully enough.™ ™ 

The second part of the criticism against QCD sum rules is more essential. It is 
concerned with the ansatz taken to parametrize the spectral function. For instance, 
it is common to assume the "pole + continuum" functional form, where the pole 
represents the hadron in question and the continuum stands for the excited and 
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scattering states. While this ansatz may be justified in cases where the low-energy 
part of the spectral function is dominated by a single pole and the continuum states 
become important only at higher energies (the /O-meson channel for instance is such 
a case), it is not at all clear if it is also valid in other cases. For example, as shown in 
Ref. [5j), where the u-meson channel was investigated using QCD sum rules, it can be 
difficult to distinguish the continuum spectrum from a broad resonance, because they 
lead to a similar behavior of the "pole" mass and residue as functions of the Borel 
mass. Moreover, using the "pole + continuum" ansatz, the outcome of the analysis 
usually depends on unphysical parameters such as the Borel mass or the threshold 
parameter, and it is not always a trivial matter to determine these parameters in a 
consistent way. After all, the ansatz of parametrizing the spectral function makes a 
full error estimation impossible in QCD sum rules. 

As a possible solution to these problems, we propose to analyze QCD sum rules 
with the help of the Maximum Entropy Method (MEM) . This method has already 
been applied to Monte-Carlo studies in both condensed matter physics^ and lattice 
QCD,^'^ and has applications in many other areas.'^l' It makes use of Bayes' theorem 
of probability theory, which helps to incorporate known properties of the spectral 
function such as positivity and asymptotic values into the analysis and finally makes 
it possible to obtain the most probable spectral function without having to introduce 
any additional a priori assumptions about its explicit form. It even allows us to 
estimate the error of the obtained spectral function. Therefore, using this approach, 
it should in principle be possible to study the spectral function of any channel, 
including those for which the "pole + continuum" assumption is not appropriate. 

However, as a first step it is indispensable to check whether QCD sum rules are 
a suitable target for MEM and if it is possible to obtain any meaningful information 
on the spectral function by this method. To provide an answer to these questions is 
the main object of this paper. To carry out this check, we have chosen to investigate 
the sum rule of the p-meson. This channel is one of the first subjects that have 
been studied in QCD sum rules and it is fair to say that it is the channel where this 
method has so far seen its most impressive success. As mentioned earlier, it is a 
case where the "pole + continuum" ansatz works well and we thus do not expect to 
gain anything really new from this analysis. Nevertheless, apart from the aspect of 
testing the applicability of our new approach, we believe that it is worth examining 
this channel once more, as MEM also provides a new viewpoint of looking at various 
aspects of this particular sum rule. 

The paper is organized as follows. In §2, the formalism of both QCD sum rules 
and MEM are briefiy recapitulated. Then, in §3, we show the findings of a detailed 
MEM analysis using mock data with realistic errors. This is followed by the results 
of the investigation of the actual sum rule of the /3-meson channel in §4. Finally, the 
summary and conclusions are given in §5. 



A Bayesian Approach to QCD Sum Rules 3 

§2. Formalism 

2.1. Basics of QCD sum rules 

The essential idea of the QCD sum rule approacbP''2J is to make use of the 
analytic properties of the two-point function of a general operator J[x): 



n{q^) = i /"d^xe^'?^(0|T[J(x)jt(0)]|0). (2-1) 



The analytic properties of this expression allow one to write down a dispersion rela- 
tion, which connects the imaginary part of n{q^) with its values in the deep euclidean 
region, where q is large and spacelike {—q^ — >• oo). This is the region, where it is 
possible to systematically carry out the operator product expansion. The dispersion 
relation reads as: 

, o, I f°° , Imn(s) , , , , , 

n(q^) = - / ds V^ + (subtraction terms). (2-2) 

T^ Jo s-q^ 

Here, the "subtraction terms" are polynomials in q^, which have to be introduced to 
cancel divergent terms in the integral on the right-hand side of Eq. (|2-2|) . When the 
Borel transformation is applied, these terms vanish and we therefore do not consider 
them any longer in the following. Using the Borel transformation, defined below, also 
has the advantage of improving the convergence of the operator product expansion. 



Lm= lim -L^i^- — . (2-3) 




After having calculated the OPE of n{q^) in the deep euclidean region, the Borel 
transformation is applied to both sides of Eq. (j2-2p , giving the following expression 
for Gope{M) = LM[n^^^{q^)\, which depends only on M, the Borel mass: 



1 /'OO 

Gope{M) = —-2 / (ise-^/*^'lm77(s) 

o /"OO 



(24) 



Here, s = w^ and Im77(s) = vr/9(a;) were used. At this point, one usually has to 
make some assumptions about the explicit form of the spectral function p{oj). The 
most popular choice is the "pole + continuum" assumption, where one introduces a 
5-function for the ground state and parametrizes the higher-energy continuum states 
in terms of the expression obtained from the OPE, multiplied by a step function: 

Imi7(s) = ^\\\^5{s - n?) + d{s - Sth)lmn^^^{s). (2-5) 

Here, |Ap is the residue of the ground-state pole and Sth determines the scale, above 
which the spectral function can be reliably described using the perturbative OPE ex- 
pression. It is important to note that Sth in many cases does not directly correspond 
to any physical threshold, which often makes its interpretation rather difficult. 
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We will in the present study follow a different path and do not use Eq. ()2-5p 
or any other parameterizations of the spectral function. Instead, we will go back to 
the more general Eq. (|2-4|) and employ the Maximum Entropy Method to directly 
extract the spectral function from this equation and the known properties of the 
spectral function at high and low energies. This procedure reduces the amount of 
assumptions that have to be made in the analysis. 

2.2. Basics of the Maximum Entropy Method (MEM) 

In this subsection, the essential steps of the MEM procedure are reviewed in 
brief. For more details, see Ref. E]), where applications in condensed matter physics 
are discussed, or Refs. [7]) and [8]) for the implementation of this method to lattice 
QCD. For the application of MEM to Dyson- Schwinger studies, see also Ref. 10]). 

The basic problem that is solved with the help of MEM is the following. Suppose 
one wants to calculate some function p{uj), but has only information about an integral 
of p{u)) multiplied by a kernel K{M,uj): 

Gp{M)= I dLoK{M,Lo)p{uj). (2-6) 

For QCD sum rules, this equation corresponds to Eq. (|2-4p with 

K(M,.) = ^e--^A-^^ (2.7) 

and GoPEiM) instead of Gp{M). When Gope{M) is known only with limited 
precision or is only calculable in a limited range of M, the problem of obtaining p{oo) 
from Gope{M) is ill-posed and will not be analytically solvable. 

The idea of the MEM approach is now to use Bayes' theorem, by which additional 
information about p{u}) such as positivity and/or its asymptotic behavior at small 
or large energies can be added to the analysis in a systematic way and by which one 
can finally deduce the most probable form of p{uj). Bayes' theorem can be written 
as 

where prior knowledge about p{uj) is denoted as H, and P[p\GH] stands for the 
conditional probability of p{uj) given GoPEiM) and H. Maximizing this functional 
will give the most probable p{uj). P[G\pH] is the "likelihood function" and P[p\H] 
the "prior probability". Ignoring the prior probability and maximizing only the 
likelihood function corresponds to the ordinary x^-fitting. The constant term P[G\H] 
in the denominator is just a normalization constant and can be dropped as it does 
not depend on p{uj). 

2.2.1. The likelihood function and the prior probability 

We will now briefly discuss the likelihood function and the prior probability one 
after the other. Considering first the likelihood function, it is assumed that the func- 
tion Gope{M) is distributed according to uncorrelated Gaussian distributions. For 
our analysis of QCD sum rules, we will numerically generate uncorrelated Gaussianly 
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distributed values for GoPEi^), which satisfy this assumption. We can therefore 
write for P[G\pH]: 



P[G\pH] =e 



'L[p] 



LM = : /'"- <iM [°°-<'"' - °''<*^« ' ''■'> 



Here, Gope{M) is obtained from the OPE of the two-point function and Gp{M) is 
defined as in Eq. ()2-6p and, hence, imphcitly depends on p{ijj). a{M) stands for the 
uncertaintity of Gope{M) at the corresponding value of the Borel mass. In practice, 
we will discretize both the integrals in Eqs. ()2-6p and (|2-9p and take A'^ data points 
for p{u}) in the range a; = ~ Wmax = 6.0 GeV and Nm data points for Gope{M) 
and Gp{M) in the range from Mmin to Mmax- 

The prior probability parametrizes the prior knowledge of p{oj) such as positivity 
and the values at the limiting points, and is given by the Shannon-Jaynes entropy 
S[p]: 



P[p\H] - 


= e"5[p]^ 


S[p\- 


/■oo 

= / du) 
Jo 



p{uj) — m{oj) — p{oj) lo~' 



m{uj) 



(2-10) 



For the derivation of this expression, see for instance Ref. [S]). It can be either derived 
from the law of large numbers or axiomatically constructed from requirements such as 
locality, coordinate invariance, system independence and scaling. For our purposes, 
it is important to note that this functional gives the most unbiased probability for 
the positive function p{uj). The scaling factor a that is newly introduced in Eq. 
(|2-lUp will be integrated out in a later step of the MEM procedure. The function 
m{uj), which is also introduced in Eq. ()2-10p . is the so-called "default model". In 
the case of no available data Gope{M), the MEM procedure will just give m{uj) for 
p{uj) because this function maximizes P[/9|i:f]. The default model is often taken to 
be a constant, but one can also use it to incorporate known information about p{uj) 
into the calculation. In our QCD sum rule analysis, we will use m{uj) to fix the value 
of p{u)) at both very low and large energies. As for the other integrals above, we will 
discretize the integral in Eq. ()2T0p and approximate it with the sum of N^^ data 
points in the actual calculation. 

2.2.2. The numerical analysis 

Assembling the results from above, we obtain the final form for the probability 

P[p\GH]: 

P[p\GH] ^P[G\pH]P[p\H] 

=e^^p\ (2-11) 

Q[p]^aS[p]-L[p]. 

It is now merely a numerical problem to obtain the form of p{u)) that maximizes 
Q[p] and, therefore, is the most probable p{uj) given Gope{M) and H. As shown 
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for instance in Ref. [8]), it can be proven that the maximum of Q[p\ is unique if it 
exists and, therefore, the problein of local minima does not occur. For the numerical 
determination of the maximum of Q[p\i we use the Bryan algorithm,'^^ which uses 
the singular-value decomposition to reduce the dimension of the configuration space 
of p{oj) and, therefore, largely shortens the necessary calculation time. We have 
also introduced some slight modifications to the algorithm, the most important one 
being that once a maximum is reached, we add to p{oj) a randomly generated small 
function and let the search start once again. If the result of the second search agrees 
with the first one within the requested accuracy, it is taken as the final result. If not, 
we add to p{uj) another randomly generated function and start the whole process 
from the beginning. We found that this convergence criterion stabilizes the algorithm 
considerably. 

Once a Pai^) that maximizes Q[p] for a fixed value of a is found, this parameter 
is integrated out by averaging p{uj) over a range of values of a and assuming that 
P[p\GH] is sharply peaked around its maximum P[pa\GH]: 



/>out(w) = / [dp] / dapioj)P[p\GH]P[a\GH] 

'' (2-12) 

dapa{uj)P[a\GH]. 



This /Oout('^) is our final result. To estimate the above integral, we once again make 
use of Bayes' theorem to obtain P[a\GH\ (for the derivation of this expression, see 
again Ref. E])): 

P[a|GF]ocP[Q|F]expr^Vlog^— + QK]]. (2-13) 

L2 ^ a + Xk J 

Here, A^ represents the eigenvalues of the matrix 



d'L 



^- = ^a^^^ 



(2-14) 



p=pa 



where pi stands for the discretized data points of p{u)): pi = p{uji)Auj, with Auj = 
^i^^^ and Wj = j^ X Wmax- To get to the final form of the right-hand side of Eq. 
(|2-13p , we made use of the fact that the measure [dp] is defined as 

In Eq. (|2-13p . we still have to specify P[a|-ff], for which one uses either Laplace's rule 
with i-*[a|i?] = const, or Jeffreys' rule with P[a|//] = — . We will employ Laplace's 
rule throughout this study, but have confirmed that the analysis using Jeffreys' rule 
gives essentially the same results, leading to a shift of the final p-meson peak of only 
10 - 20 MeV. 

To give an idea of the behavior of the probability P[a[GH], a typical example of 
this function is shown in Fig. [TJ Its qualitative structure is the same for all the cases 
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Fig. 1. Probability P[a\GH], as obtained in our analysis of mock data. To generate this function, 
the errors of Ref. 15.) and the default model of Eq. p-3|) with parameters uio — 2.0 GeV and 
S — 0.1 GeV have been used. The dashed line corresponds to the value of 0.1 x PfamaxIG-ff], 
from which the boundaries of the integration region, qiow and Qup, are determined. 



studied in this paper. To calculate the integral of Eq. (|2-12p . we first determine the 
maximum value of P[a\GH], which has a pronounced peak: P[aaiax\GH]. Then, we 
obtain the lower and upper boundaries of the integration region (aiow, ceup) from the 
condition P[a\GH] > 0.1 x P[ama.x\GH] and normalize P[a\GH] so that its integral 
within the integration region gives 1. After these preparations, the average of Eq. 
()2-12p is evaluated numerically. 

2.2.3. Error estimation 

As a last step of the MEM analysis, we have to estimate the error of the obtained 
result /?out('^)- The error is calculated for averaged values of po^t{^) over a certain 
interval (wi, L02), as shown below. 

The variance of p{uj) from its most probable form for fixed a, Spaiuj) = p(^) — 

Pa{<^)i averaged over the interval {ijJi^ ijJ2)-, is defined as 

1 r r^2 

[dp] I dujduj'5pci{oj)6pa{u:')P[p\GH] 



{{5po.y 



'Wl,Ld2 



(W2 



■UJi) 
1 



(a;2 -ujiY 



W2 



dujduo' 



wi 



5^Q 



5p{ijj)5p{ijj'] 



(2-16) 



P=Pc. 



where again the definition of Eq. (|2-15p and the Gaussian approximation for P[p\GH] 
were used. Taking the root of this expression and averaging over a in the same way 
as was done for Paiyj)-, we obtain the final result of the error 5pout(w), averaged over 
the interval (cji, a;2): 



{5P. 



OVit/LjJ\,U2 



"lo 



da^J{{6po 



,P[a\GH]. 



(2-17) 



The interval (wi, UJ2) is usually taken to cover a peak or some other structure of 
interest. The formulas of this subsection will be used to generate the error bars of 
the various plots of /Ooutl"^)) as shown in the following sections. 
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§3. Analysis using mock data 

The uncertainties that are involved in QCD sum rule calculations mainly origi- 
nate from the ambiguities of the condensates and other parameters such as the strong 
coupling constant or the quark masses. These uncertainties usually lead to results 
with relative errors of about 20%. It is therefore not a trivial question if MEM can 
be used to analyze the QCD sum rule results, or if the involved uncertainties are too 
large to allow a sufficiently accurate application of the MEM procedure. 

To investigate this question in detail, we carry out the MEM analysis using mock 
data and realistic errors. Furthermore, we will study the dependence of the results 
on various choices of the default model m{(jj) and determine which one is the most 
suitable for our purposes. This analysis will also provide us with an estimate of the 
precision of the final results that can be achieved by this method and what kind of 
general structures of the spectral function can or cannot be reproduced by the MEM 
procedure. 

3.1. Generating mock data and the corresponding errors 

Following Refs. [S]) and ll2|) . we employ a relativistic Breit-Wigner peak and a 
smooth function describing the transition to the asymptotic value at high energies 
for our model spectral function of the p-meson channel: 

( ^ _ ^^P ^P^P I ^ /"l I "^^ ^ 

^'^^"^^ vr (a;2-m2)2+r2m2 ^47r2V ^ ^;i + e(^o— ')/5' 

The values used for the various parameters are 

nip = 0.77 GeV, m^ = 0.14 GeV, 
wo = 1-3 GeV, 6 = 0.2 GeV, 
gpn^ = 5.45, as = 0.5, (3-2) 

Fp = ^^ = 0.141 GeV. 

9p7V7T 

The spectral function of Eq. (j3-l|) is then substituted into Eq. (|2-4|) and the integra- 
tion over cj is performed numerically to obtain the central values of the data points 
of Gmock(-^)- (In this section, we will use Gynock{M) instead of Gope{M) to make 
it clear that we are analyzing mock data.) 

We now also have to put some errors a{AI) to the function Gmock(-^)- To make 
the analysis as realistic as possible, we will use exactly the same errors as in the actual 
investigation of the OPE results. How these errors are obtained will be discussed 
later, in the section where the real OPE results are analyzed. We just mention here 
that when analyzing the OPE results, we will use three different parametrizations 
for the condensates and other parameters, namely, those given in Refs. [T3 ]) - I15p (see 
Table |I|). These parametrizations lead to different estimations of the errors, but 
for the mock data analysis of this section, these differences are not very important. 
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Fig. 2. Central values of the mock data obtained by numerically integrating Eq. (|2-4|) with Eq. 
p-l|) . The errors of the data, extracted from the parameters of Ref. ,15;), are indicated by the 
gray region. The lower boundary of the shown errors corresponds to Afmin, below which the 
OPE does not converge. For comparison, the integral of Eq. H2-4|) is also shown for the case 
when only the continuum part, the second term of Eq. p-ip . is taken for the spectral function. 



Here, we will therefore mainly use the errors obtained from the parameters of Ref. 
lisp . The resulting function Gmock(-^) is given in Fig. [21 together with the range 
G„,ock(M)±a(M). 

3.1.1. Determination of the analyzed Borel mass region 

Next, we have to decide what range of M to use for the analysis. For the lower 
boundary, we can use the usual convergence criterion of the OPE such that the 
contribution of the highest-dimensional operators is less than 10% of the sum of all 
OPE terms. This is a reasonable choice, as the errors originating from the ranges 
of condensate values lead to uncertainties of up to 20%, and it would therefore not 
make much sense to set up a more strict convergence criterion. For the parameters 
of Ref. IS]), this criterion gives Mmin = 0.77 GeV. 

Considering the upper boundary of M, the situation is less clear. In the conven- 
tional QCD sum rule analysis, it is standard to use the pole dominance condition, 
which makes sure that the contribution of the continuum states does not become 
too large. As we do not resort to the "pole -|- continuum" ansatz in the current 
approach, the pole dominance criterion does not have to be used and one can, in 
principle, choose any value for the upper boundary of M. Nevertheless, because we 
are mainly interested in the lowest resonance peak, we will use a similar pole dom- 
inance criterion as in the traditional QCD sum rules. By examining the mock data 
in Fig. [21 one sees that while the resonance pole contributes most strongly to the 
data around M ~ 1 GeV, the contributions from the continuum states grow with 
increasing M and finally start to dominate the data for values that are larger than 
1.5 GeV. We will therefore use Mmax = 1-5 GeV as the upper boundary of M for 
the rest of the present paper. The dependence of the final results on this choice is 
small, as will be shown later in Fig. [5l 

Finally, the values of the Nm data points of G^ocki^) between Mmin and 
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Mmax are randomly generated, using Gaussian distributions with standard devia- 
tions a{M), cenered at the values obtained from the integration of Eq. (j2-4p . The 
ranges of values of G^oc\{M) are indicated by the gray region in Fig. [21 We take 100 
data points for functions of the Borel mass M {Nm = 100) and have checked that 
the results of the analysis do not change when this value is altered. For functions of 
the energy w, 600 data points are taken {N^ = 600). 

3.2. Choice of an appropriate default model 

It is important to understand the meaning of the default model m(a;) in the 
present calculation. It is in fact used to fix the value of the spectral function at high 
and low energies, because the function GoPEiM) contains only little information on 
these regions. This can be understood firstly by considering the property of the kernel 
of Eq. (|2-7|) . which is zero at w = 0. Gope{M) therefore contains no information 
on pin{uj = 0) and the corresponding result of the MEM analysis Poutl"^ = 0) will 
thus always approach the default model, as Gope{M) does not constrain its value. 
Secondly, we use Gope{M) only in a limited range of M, because the operator 
product expansion diverges for small M and the region of very high Gope{M) is 
dominated by the high-energy continuum states, which we are not interested in. 
The region of the spectral function, which contributes most strongly to Gope{M) 
between M^i^ and Mmax; lies roughly in the range between Wmin (— -^min) and Wmax 
(~ Mjnax)i as can be for instance inferred from Fig. [51 pout{^) will then approach 
the default model quite quickly outside of this region, because there is no strong 
constraint from Gope{^)- This implies that the values of /Oout(w) at the boundaries 
are fixed by the choice of the default model and one should therefore consider these 
boundary conditions as inputs of the present analysis. Once these limiting values 
of p{oj) are chosen, the MEM procedure then extracts the most probable spectral 
function po\it{<^) given Gope{M) and the boundary conditions supplied by m{ijj). 

To illustrate the importance of choosing appropriate boundary conditions, we 
show the results of the MEM analysis for a constant default model, with a value fixed 
to the perturbative result at high energy. Here, the boundary condition for the low 
energy is not correctly chosen, because the spectral function is expected to vanish 
at very low energy. The result is given in Fig. [3] and clearly shows that Pout(^) does 
not reproduce p\n{uj). 

This is in contrast to the corresponding behavior in lattice QCD, where it suffices 
to take a constant value of the spectral function, chosen to be consistent with the 
high-energy behavior of the spectral function to obtain correct results. The reason 
for this difference is mainly that the OPE in QCD sum rules is not sensitive to the 
low-energy part of the spectral function, owing to the properties of the kernel and our 
limitations in the applicability of the OPE. Most importantly, the information that 
there is essentially no strength in the spectral function below the rho-meson peak, is 
stored in the region of the Borel mass M around and below M = 0.5 GeV. However, 
as the OPE does not converge in this region, it is not available for our analysis and 
we therefore have to use the default model to adjust the spectral function to the 
correct behavior. On the other hand, in lattice QCD, it is possible to calculate the 
correlator reliably at sufficiently high euclidean time, where the effective mass plot 
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Fig. 3. Outcome of the MEM analysis using a constant default model with its value fixed to the 
perturbative result, pin(i^) is the function that was used to produce the mock data, given in Eq. 
(|3-ip . and poMti'-'j) shows the result of the MEM procedure. 



reaches a plateau and thus only the ground state contributes. Therefore, it seems 
that from lattice QCD one can gain sufficient information on the low-energy part of 
the spectral function, and one does not need to adjust the default model to obtain 
physically reasonable results. 

For the present analysis, we introduce the following functional form, to smoothly 
connect low- and high-energy parts of the default model. 



m{uj) 



47r2 V^ vr ; 1 



1 



(3-3) 



which is close to at low energy and approaches the perturbative value 1 + — at high 
energy, changing most significantly in the region between ujq — 5 and ojq + 6. This 
function can be considered to be the counterpart of the "continuum" in the "pole 
-|- continuum" assumption of Eq. (|2-5|) . where 5 is essentially taken to be and 
the threshold parameter sth corresponds to ujq. Equation ()3-3p nevertheless enters 
into the calculation in a very different way than the second term of Eq. (j2-5p in the 
conventional sum rules, and one should therefore not regard these two approaches 
to the continuum as completely equivalent. 

We have tested the MEM analysis of the mock data for several values of ojq 
and 5 and the results are shown in Fig. HI One sees that in the cases of d) and 
e), the sharply rising default model induces an artificial peak in the region of uoq. 
Even though these peaks are statistically not significant, they may lead to erroneous 
conclusions. We will therefore adopt only default models, for which only small 
artificial structures appear, such as in the case of the figures a), b) and c). Comparing 
these three figures, it is observed that the error of the spectral function relative to 
the height of the p- meson peak is smallest for the parameters of c). We therefore 
adopt this default model with ujq = 2.0 GeV and (5 = 0.1 GeV in the following 
investigations. 

It is worth considering the results of Fig. [J] also from the viewpoint of the 
dependence of the peak position on the default model. It is observed that even 
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Fig. 4. Results of the MEM investigation of mock data with various default models. As in Fig. [3l 
the solid lines stand for the output of the analysis /9out(i^), the long-dashed lines for the default 
model with the parameters shown in the figure, and the short-dashed lines for the input function 
pin(i^) of Eq. p-l|l . The horizontal bars show the values of the spectral function, averaged over 
the peaks {pout)ui,uj2 a-^d the corresponding ranges {pout)u,i,LJ2 ± (5 pout) 1^1,1^2- Their extent 
shows the averaged interval (a;i, 1^2). For figures c), d) and e), the lower error bars of the second 
peak are not shown because they lie below p(a;) = 0. 



though the height of the p-meson peak varies quite strongly depending on which 
defauh model is chosen, its position only varies in a range of ±40 MeV, which shows 
that the present evaluation of the lowest pole position does not strongly depend on 
the detailed values of ujq and 6. This behavior should be compared with the results 
of the usual sum rules, where the dependence of the pole mass on the threshold 
parameter Sth is stronger. On the other hand, we have to mention that wq should 
not be chosen to have a value much below ujq = 1.5 GeV, because, in this region, 
the artificial structures discussed above start to interfere with the /o-meson peak. 
Moreover, if the default model approaches the limit shown in Fig. [3l where m(a;) is 
just a constant fixed to the asymptotic value at high energy, the /j-meson peak will 
gradually disappear. 

3.3. Investigation of the stability of the obtained spectral function 

Next, we will briefly discuss the dependence of our results on the upper boundary 
of the employed Borel mass region Mmax- In Fig- El we show the results for the values 



M„ 



1.5 GeV, 2.0 GeV and 2.5 GeV. Here, the default model with parameters 
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bjQ = 2.0 GeV and (5 = 0.1 GeV was used. The spectral functions of these three 
cases almost coincide and have the same qualitative features. Quantitatively, the 
peak position of the p-meson is shifted only 20 MeV when Mmax is raised from 1.5 
GeV to 2.5 GeV. 

Let us also check how the results of the analysis are affected by a different 
choice of parameters, leading to altered magnitudes of error and also differing lower 
boundaries of the Borel mass Mmin- Of the three parameter sets used in this study, 
given in Table [H the errors of Ref. I13p are rather small, while the errors of Refs. 
[T^ and[T5|) are larger and have about the same overall magnitude. Moreover, Mminj 
which is determined from the OPE convergence criterion mentioned earlier, takes 
values Mmin = 0.72 GeV for Ref. [13]), 0.83 GeV for Ref. [H]) and 0.77 GeV for Ref. 
[T5]) . To understand how these parameters affect the MEM analysis, the results of 
the calculation using the same central values for the mock data, but different errors 
and Mmiii) are shown in Fig. [6l It is observed that the spectral functions are very 
similar and depend only weakly on the choice of errors and Mmin. 

Finally, it is important to confirm whether the lowest peak that we have observed 
in all the results so far really originates from the p-meson peak of the input spectral 
function. In other words, we have to verify that the lowest peak obtained from the 
MEM analysis really disappears when the p-meson peak is removed from the input 
spectral function. The result for this case is given in Fig. [71 We see that while 
we get the same (nonsignificant) peak around 2.0 GeV as before, which is induced 
by the sharply rising default model in this region, the lower peak has completely 
disappeared. This confirms that the lower peak is directly related to the /9-meson 
peak and is not generated by any other effects of the calculation. 

3.4. Estimation of the precision of the final results 

To obtain an estimate of the precision of the current approach, let us now turn 
to the numerical results of our analysis of mock data. We regard plot c) of Fig. 
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Fig. 5. MEM results for three different values of the upper boundary of the Borel mass Mmax- For 
the default model, the same version as in part c) of Fig. [4]was used. Thus, the solid line of this 
figure is the same as the one of part c) in Fig. U 
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Fig. 6. Findings of the MEM analysis with different error estimations and lower boundaries of the 
Borel mass i\fmin- For the default model, the same version as in part c) of Fig. [4] was used. 
The solid line uses parameters of Ref. 1131 , the long-dashed lines those of Ref. I14p and the 
short-dashed line those of Ref. I15|l . Note that we use for this plot the same mock data for all 
three cases and vary only the errors and Mmin- 
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Fig. 7. Outcome of the MEM analysis using mock data without the p-meson peak. pin{(jj) is the 
function that was used to produce the mock data and pout(i^) shows the result of the MEM 
procedure. For the default model, the same version as in part c) of Fig. 2] was used. 



m as our main result, the central value of the peak being nip^out = 0.84 GeV. The 
shift from the true value of mp^[^ = 0.77 GeV is caused by the errors of the involved 
parameters and by the fact that we cannot use all the data points of G^ocki^), but 
only the ones for which we have some confidence that the OPE converges. 

Furthermore, as discussed above, there are some additional uncertainties of ±40 
MeV coming from the choice of the default model and ±20 MeV from the value of 
-^max- The overall error is then obtained by taking the root of the sum of all squared 
errors and rounding it up. This gives 

Z\mp~ ±90 MeV, (3-4) 

which is quite large but seems to be realistic when one considers the large errors of 
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Fig. 8. An example of a result of the fitting procedure described in the text. For the spectral 
function pout(^), we use here the same as in part c) of Fig. |4l For the peak, a relativistic Breit- 
Wigner form is employed, while the background is parametrized by a second-order polynomial 
with positive valued coefficients. 



the condensates that are involved in the calculation. For the other parametrizations, 
we get similar results, concretely, Arup ~ ±60 MeV for Ref. IT3]l and Arup ~ ±100 



MeV for Ref. [E 

Having the spectral function at our disposal, it becomes possible to extract 
the coupling strength of the interpolating field to the p-meson state, the quantity 
corresponding to Fp in Eq. ()3-ip . We obtain this coupling strength by fitting the 
spectral function in the region of the /o-meson resonance with a relativistic Breit- 
Wigner peak of the same functional form as the first term of Eq. (|3-ip plus a 
second-order polynomial to describe the continuum background. In order that the 
background does not become negative and does not contribute to the peak, we have 
restricted the coefficients of the polynomial to positive values. An example of the 
resulting curves is given in Fig. [5J For the case of plot c) in Fig. |H we get a value 
of Fp^ont = 0.178 GeV, which is somewhat larger than the true value of -Fp,in = 0.141 
GeV. It is not a surprise that the precision of this quantity is worse than that of 
the peak position, because the shape of the peak is deformed quite strongly owing 
to the MEM procedure. As can be suspected when looking at Fig. UJ there is also a 
fairly large uncertainty coming from the choice of the default model, which is about 
±0.031 GeV. On the other hand, we found that the dependence on the boundaries 
of the fitting region and on Mmax is very small, being ±0.003 GeV and ±0.001 GeV, 
respectively. Altogether, this gives the following error for Fp: 



Z\F„~ ±0.049 GeV. 



(3-5) 



A similar analysis for the parameters of Ref. I13p gives AFp ~ ±0.038 GeV and 
AFp ~ ±0.049 GeV for Ref. [HI). 

As one further point, it is important to investigate if and how the precision of the 
MEM analysis improves once the OPE data will be available with better precision. 
To answer this question, we have repeated the analysis using an error with a smaller 
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Fig. 9. Two versions of G'mock(Af), the solid line showing the one with the standard value for the 
width as given in Eqs. (|3-ip and (I3-2[I . and the dashed line giving a version where Fpiuj) of Eq. 
(|M|) is doubled and Fp = 0.149 GeV is used. 



magnitude and have found that, concerning the pole position, the reproducibihty 
indeed improves with a smaller error. For instance, using the errors obtained from 
the loffe parameters of Ref. [T5]l . we get mp^out = 0.84 GeV, but when we reduce 



the errors by hand to 



of their real value, the result shifts to m 



p,out 



0.78 



GeV, which almost coincides with the correct value of mp^n = 0.77 GeV. To a lesser 
degree, the same is true for the residue. Its value changes from Fp^out = 0.178 GeV 
to -Fp,out = 0.167 GeV with reduced error values, compared with the input value 
of -Fp,in = 0.141 GeV, which is also an improvement. On the other hand, we could 
not observe a significantly better reproducibility of the width with the reduced error 
values. 

These results show that the MEM analysis of QCD sum rules has the potential 
to become more accurate in obtaining the position and the residue of the p-meson 
pole, once the values of the condensates are known with better precision. At the 
same time, it has to be noted that an accurate determination of the width seems to 
be difficult to achieve in the current approach even with smaller errors. The reason 
for this is discussed in the following section. 

3.5. Why it is difficult to accurately determine the width of the p-meson 

We have so far focused our discussion on the reproduction of the parameters 
nip and Fp of the spectral function and have shown that, by the MEM analysis, 
they could be reproduced with a precision of approximately 20%. Considering now 
the width Fp, the situation turns out to be quite different, as can be observed for 
instance in Fig. HI We see there that the values of the width come out about twice 
as large as in the input spectral function, Eq. (|3Tp . 

The reason for this difficulty of reproducing the width lies in the small depen- 
dence of Gp{M) on Fp. This is illustrated in Fig. [9l We see there that the curve 
obtained from our model spectral function of Eqs. (|3-ip and (|3-2p and the one gen- 
erated from a spectral function for which the width of the /?-meson peak has been 
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doubled and a slightly larger value for Fp has been used {Fp = 0.149 GeV compared 
with the standard value of 0.141 GeV), almost coincide. This means that with the 
precision available for our QCD sum rule analysis, it is practically impossible to 
distinguish between these two cases. Examining the curves of Fig. [9] a bit more 
carefully, it is found that the most prominent difference between them lies in the 
region of 0.5 GeV or below. However, this region cannot be accessed by the current 
calculation, as the OPE does not converge well for such small values of M. Further- 
more, even if we could have calculated the OPE to higher orders and would thus 
have some knowledge about Gope{M) in the region below 0.5 GeV, this would most 
likely not help much, as the uncertainty here will be large owing to the large number 
of unknown condensates that will appear at higher orders of the OPE. We therefore 
have to conclude that it is not possible to say anything meaningful about the width 
of the /9-meson peak at the current stage. To predict the width, the OPE has to be 
computed to higher orders and the various condensates have to be known with much 
better precision than they are today. 

§4. Analysis using the OPE results 

4.1. The p-meson sum rule 

Carrying out the OPE and applying the Borel transformation, we obtain Gope{M), 
the left-hand side of Eq. ()2-4p . In the case of the vector meson channel, we use the 
operator 

i^i = i^{ui^iU-d-i^d), (4-1) 



which stands for J in Eq. ()2-ip and take the terms proportional to the structure 
QiJ-Qy ~ Q^9nu- We then arrive at the following expression for GoPEiM), where the 
OPE has been calculated up to dimension 6: 

GoPEiM) ~(l + ^{as)) + (2m{qq) + ^{^G'))^^ 



""asKiqqf— : + ..., (4-2) 



81 ^ ^^^' M6 
rj{as) =— + 0.154a2 _ 0.372a^ + . . . 



vr 



Here, as is the usual strong coupling constant, given as |^, ?tt, stands for the (av- 
eraged) quark mass of the u- and d-quarks, and (qq) is the corresponding quark 
condensate. Meanwhile, the gluon condensate l^G'^) is an abbreviated expression 
for ^^G" G"'^^) and k parametrizes the breaking of the vacuum saturation approx- 
imation, which has been used to obtain the above result for n = 1. 

A few comments about this result are in order here. For the perturbative term, 
which is known up to the third order in as, we have taken the number of flavors to 
be Nf = 4. Note that only the second and third terms of rj(as) depend (weakly) on 
Nf^ and that the final results of the analysis are thus not affected by this choice. 
We have not considered the running of ag in deriving Eq. (j4-2p for simplicity. If 
the running is taken into account, the coefficient of the third term of rj^ag) changes 
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due to the Borel transformation, as was shown in Ref. [T7|l . Nevertheless, this again 
leads only to a minor change in the whole expression of Eq. (j4-2p and does not alter 
any of the results shown in the following. Considering the terms proportional to 
1/M^, the first-order corrections of the Wilson coefficients are in fact known,'^ but 
we have not included them here as the values of the condensates themselves have 
large uncertainties, and it is therefore not necessary to determine the corresponding 
coefficients with such a high precision. It is nonetheless important to note that these 
corrections are small (namely of the order of 10% or smaller, compared with the 
leading terms) and thus do not introduce any drastic changes into the sum rules. 

4.1.1. Values of the parameters and their uncertainties 

There are various estimates of the values of the condensates and their ranges. 
We will employ the ones given in three recent publications: [T3 l) -tl5 p . The explicit 
values are given in Table HI where they have been adjusted to the renormalization 
scale of 1 GeV. For the central value of m{qq), we make use of the Gell-Mann-Oakes- 
Renner relation, which gives m{qq) = —■^m'^f^ and take the experimental values of 
TTi-TT and /tt for all three cases, leading to 

m{qq) = -8.5 x 10"^ GeV^. (4-3) 



Note that due to its smallness, the term containing m{qq) does not play an important 
role in the sum rules of Eq. ()4-2p . The values of Table |T] agree well for {qq) and a^, 
while they differ considerably for (^G'^\ and n. Namely, Ref. [T^ employs values 
for (^G'^) and k that are about two times larger than those of Refs. IT3|) and 
[T5|) . Considering the error estimates of the parameters, Ref. [T3|) uses altogether 
the smallest errors as the breaking of the vacuum saturation approximation is not 
considered and only a fixed value for Ug is employed. Comparing the results obtained 
from these three parameter sets will provide us with an estimate of the order of the 
systematic error inherent in the current calculation. 

4.1.2. Determination of the errors of Gope{M) 

As can be inferred from Eq. (j4-2p and Table [H the uncertainty of Gope{M) 
will vary as a function of M and will be larger for small values of M because the 
contributions of the higher-order terms with large uncertainties of the condensates 
become more important in that region. To accurately estimate this error, we follow 
Ref. [T8|l and numerically generate Gaussianly distributed values for the condensates, 
and then examine the distribution of the resulting values of Gope{^)- For illustra- 
tion, we show the values of Gope{M = 0.84 GeV) in Fig. [TOl where the parameter 



Table I. Values of the parameters used in the calculation. These have been adjusted to the renor- 
malization scale of 1 GeV. 

Colangelo and Khodjamirian^^' NarisorP^ lofffU^ 

{qq) -(0.24 ±0.01)^ GeV^ -(0.238 ± 0.014)^ GeV^ -(0.248 ± 0.013)^ GeV^ 

(^G^) 0.012 ± 0.0036 GeV** 0.0226 ± 0.0029 GeV 0.009 ± 0.007 GeV* 

k" 1 2.5 ±0.5 1.0 ±0.1 

cts 0.5 0.50 ±0.07 0.57 ±0.08 
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estimates of Ref. lb) have been used. (t{M), the error of Gope{M), can be easily 
obtained from the formula of the standard deviation of a given set. 

As for the analysis of the mock data, the data points of Gope{M) are randomly 
generated, using Gaussian distributions with standard deviations o"(M), centered at 
Gope{M) of Eq. ([¥2]). We here again take Nm = 100 and N^ = 600. Mmin is 
determined from the 10% convergence criterion and Mmax is fixed at 1.5 GeV. 

4.2. Results of the MEM analysis 

Having finished all the necessary preparations, we now proceed to the actual 
MEM analysis of the real OPE data. First, we show the central values of the right- 
hand side of Eq. (|4-2p and the corresponding errors for the three parameter sets 
of Refs. [T3|)4T5|) on the left side of Fig. [TTJ Comparing these figures with Fig. [2l 
we see that the OPE results and the mock data obtained from Eq. (j3-ip exhibit a 
very similar behavior, even in the region smaller than Mmin, below which we have 
no control over the convergence of the OPE. 

Using these data, we have carried out the MEM analysis. For the default model, 
we have adapted the function 



m{ijj) 



i^(l + ^(«^ 



1 



1 _|_ gCti^O-w)/!? 



(4-4) 



with parameters ojq = 2.0 GeV and 6 = 0.1 GeV, which we have found in our 
investigation of the mock data to give the best reproduction of the p-meson peak 
with the smallest relative error and no large artificial peaks. The results are shown 
on the right side of Fig. [TTJ We clearly see that all three data sets give a significant 
lowest peak, which corresponds to the /j-meson resonance. To determine the peak 
position, we simply adopt the value, where the peak reaches its highest point. The 
uncertainty of this quantity has already been estimated in our mock-data analysis 
and we employ the number that has been obtained there to specify the error of our 
final results that are given in the first line of Table HH 



16000 



3 
c 

3 
O 

o 




0.026 0.027 0.028 0.029 0.03 0.031 0.032 
Gope(M=0.84 GeV) 



Fig. 10. An example of the distributed values of Gopb(M) for M = 0.84 GeV. Here, the parameters 
of Ref. 15) have been used. 
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Next, fitting tlie spectral functions of Fig. [TT] to a relativistic Breit-Wigner 
peak plus a second-order polynomial background, as was done in Fig. [8l we have 
determined the coupling strength Fp from our obtained spectral function, leading 
to the results given in the second line of Table [Til As could be expected from our 
experience with the mock data, we get results that are all somewhat larger than the 
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Fig. 11. On the left, the central values of Gope{M) with the errors a{AI) are given. The lower 
boundary of the shown errors corresponds to Mmin, below which the OPE does not converge. 
These plots should be compared with Fig. [J] On the right, the results of the MEM analysis using 
these OPE data are displayed. The dashed lines show the default inodel, and the horizontal 
bars stand for the values of the spectral function, averaged over the peaks {pout}uii,ui2 ^''^^ 
the corresponding ranges {pout)uii,ui2 ^ {Sf3out}uii,ui2- Their extent shows the averaged interval 
{uJi,iJ2)- For the lower two figures the lower error bars of the second peak are not shown because 
they lie below p{ij) = 0. 
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Table II. Final results for the three parameter sets. The corresponding errors were determined 
from our analysis of mock data in the previous section. 
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loffe' 



Experiment 



0.75 ± 0.06 GeV 
0.172 ± 0.038 GeV 



0.84 ±0.10 GeV 
0.190 ± 0.049 GeV 



0.83 ± 0.09 GeV 
0.186 ±0.049 GeV 



0.77 GeV 
0.141 GeV 
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Fig. 12. In the left figure, the mass nip is shown as a function of the gluon condensate (^^G^). 
Here, the value of {qq) = —(0.24)^ GeV'^ has been used for the quark condensate. In the right 
figure, rap is shown as a function of the quark condensate {qq) for three different values of the 
gluon condensate. To obtain these plots, the errors of Ref. HSj have been used. 



experimental value. 

4.2.1. Dependence of the p-meson peak on the values of the condensates 

Looking at Tables |T] and [III it is interesting to observe that even though the 
parameter sets of Refs. I14|) andllSp are quite different, they lead to almost identical 
results. This fact can be explained from the dependences of the properties of the 
/9-meson resonance on {qq), (^G^) and k, as will be shown in this subsection. 

Investigating the relation between the p-meson resonance and the condensates 
is also interesting in view of the behavior of this hadron at finite temperature or 
density, as the values of the condensates will change in such environments. This 
will in turn alter the QCD sum rule predictions for the various hadron properties. 
A detailed study of this kind of behavior of the p-meson is nevertheless beyond 
the scope of the present paper and is left for future investigations. Here, we only 
discuss the change in the mass of the /o-meson when the values of the condensates 
are modified by hand. The behavior of the peak position rup is shown in Fig. [T2j For 
obtaining these results, we have used the errors of Ref. [THj) and the corresponding 
values of k and as, but have confirmed that the qualitative features of Fig. [12] 
do not depend on the explicit values of those parameters. It is seen that while 
nip decreases somewhat when (— G^) increases, its value grows quite strongly when 
[qq) increases, irrespective of the value of the gluon condensate. We found that the 
coupling strength Fp shows the same qualitative behavior, slightly decreasing with 
increasing ( — G^), and strongly increasing with increasing {qq). A similar tendency 
for nip and Fp has also been observed in Ref. ilSp . This result shows that the quark 
condensate plays an essential role in determining the properties of the p-meson. 
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It is important to note here that the correlation between nip and {qq) to a large 
part occurs due to the last term in Gope{M) of Eq. (j4-2p . which contains the 
squared quark condensate. This means that a similar (but weaker) correlation exists 
between nip and k, which is also present in the last term of Gope{M). Hence, we can 
now understand why the parameters of Refs. fHj) and llSp give such similar results: 
while the large value of the gluon condensate of Ref. [T^ should lead to a smaller nip, 
this effect is compensated by the large value of k, which shifts the mass upwards. 
Therefore, the sum of these changes cancel each other out to a large degree, the net 
effect being almost no change in the spectral function for both cases. 

§5. Summary and conclusion 

We have applied the MEM technique to the problem of analyzing QCD sum 
rules. Using MEM has the advantage that we are not forced to introduce any explicit 
functional form of the spectral function, such as the "pole + continuum" ansatz that 
has often been employed in QCD sum rule studies. This therefore allows us to 
investigate any spectral function without prejudice to its actual form. Furthermore, 
with this technique, we have direct access to the spectral function without the need 
for interpreting quantities that depend on unphysical parameters such as the Borel 
mass M and the threshold parameter sth- 

To check whether it is really possible to apply MEM to QCD sum rules, we 
have investigated the vector meson channel in detail, first with mock data obtained 
from a realistic model spectral function, and then with the actual Borel-transformed 
results of the operator product expansion. The main results of this investigation are 
as follows: 

• Most importantly, demonstrating that it is possible to extract a significant peak 
in the spectral function, which corresponds to the p-meson resonance, we could 
show that the MEM technique is quite useful for analyzing QCD sum rules. For 
both mock or OPE data, we were able to reproduce the experimental p-meson 
mass nip with a precision of about 10% and the coupling strength Fp with a 
precision of about 30%. 

• We have found that owing to the properties of the kernel of Eq. ()2-7p appearing 
in QCD sum rules, the default model ni{oj) has to be chosen according to the 
correct behavior of the spectral function at low energy. We therefore have taken 
a default model that tends to zero at w = 0. On the other hand, to give the 
correct behavior at large energies, the same default model is constructed to 
approach the perturbative value in the high-energy region. 

• The position of the p-meson peak in the spectral function has turned out to be 
quite stable against changes in various parameters of the calculation, such as 
the details of the default model or the range of the analyzed Borel mass region. 
We have shown that changing these parameters leads to a fluctuation of the 
peak position of only 20 - 50 MeV. 

• Concerning the width of the lowest lying peak, we are unable to reproduce 
the value of the input spectral function of the mock data with any reasonable 
precision and have shown that the reason for this difficulty comes from the 
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insufficient sensitivity of the data Gp{M) on the detailed form of the p- meson 
peak. To accurately estimate the width, one needs not only to go to higher 
orders in the OPE, but also has to have much more precise information on the 
values of the condensates than is available today. 
• Accompanied by a steep rise in the default model m{oj), we have observed 
the appearance of artificial peaks in the output spectral function of the MEM 
analysis. These peaks are MEM artifacts and one has to be careful not to 
confuse them with the actual peaks that are predicted using the data GoPEiM). 
These results are promising and encourage us to apply this approach to other chan- 
nels, including baryonic ones. It would also be interesting to apply this method to 
the behavior of hadrons at finite density or temperature, as it would become possible 
to directly observe the change in the spectral function in hot or dense environments 
by this approach. Furthermore, it would be of interest to study the various ex- 
otic channels, containing more than three quarks. In these channels, the scattering 
states presumably play an important role and the bayesian approach of this paper 
could help clarify the situation and provide a natural way of distinguishing genuine 
resonances from mere scattering states. 

Even though these are interesting subjects for future studies, we want to em- 
phasize here that the uncertainties involved for each channel can differ considerably 
and one thus should always carry out a detailed analysis with mock data for each 
case, before investigating the actual sum rules. This procedure is necessary to check 
whether it is possible to obtain meaningful results from the MEM analysis of QCD 
sum rules. 
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